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Abstract. This paper gives a survey of wave dynamics in the Kerr space- 
time geometry, the mathematical model of a rotating black hole in equilib- 
rium. After a brief introduction to the Kerr metric, we review the separability 
properties of linear wave equations for fields of general spin s = 0, ^, Ij 2, cor- 
responding to scalar, Dirac, electromagnetic fields and linearized gravitational 
waves. We give results on the long-time dynamics of Dirac and scalar waves, 
including decay rates for massive Dirac fields. For scalar waves, we give a 
rigorous treatment of superradiance and describe rigorously a mechanism of 
energy extraction from a rotating black hole. Finally, we discuss the open 
problem of linear stability of the Kerr metric and present partial results. 
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1. Introduction to General Relativity and Black Holes 

It is apparently the Reverend John Michell who was the first to suggest in a 
paper communicated to the Royal Society of London in November 1783 that if one 
assumed that light consisted of particles subject to the law of universal gravitation, 
then the Universe could in principle contain stars that would be completely invisible 
to any external observer. This conclusion was arrived at independently by Laplace 
in 1796 in his Exposition du systeme du monde, and the idea is credited to him in 
several classical references in General Relativity, such as The Large Scale Structure 
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of Space- Time by Hawking and Ellis [33jj. The concept of a dark star would have 
probably remained as a mathematical curiosity had it not been for the advent of the 
Theory of Relativity and for the important advances in the study of stellar struc- 
ture which took place in the early decades of the twentieth century, following the 
development of quantum statistical mechanics. These ideas led a nineteen year old 
Indian student, while on his way in July of 1930 from Madras, India, to Cambridge, 
England, where he was to begin his doctoral studies, to discover the extraordinary 
consequence of General Relativity that a star which is sufficiently massive will, after 
having exhausted its nuclear fuel, collapse gravitationally to a point at which the 
very structure of space-time will become singular. In spite of the severe opposition 
expressed to this daring work by Sir Arthur Eddington, the leading astronomer of 
his time, the Indian student's work was completely validated by further analysis 
and he went on to receive a Nobel Prize many years later, in 1983, for his discovery 
of this lower mass limit. His name was Subrahmanyan Chandrasekhar, and he in 
turn became one of the most eminent theoretical astrophysicists of his time. 

It took a long period of maturation before a proper understanding was developed 
of what the space-time geometry would be near such a collapsed object. However, it 
was already apparent from the early work dating back to 1915 of Karl Schwarzschild 
on the relativistic space-time geometry near a stationary spherically symmetric star, 
that the singularity should not be directly visible to an external observer who would 
be at rest with respect to the star, and that it should be hidden behind a a one-way 
membrane, known today as an event horizon. This led the physicist John Wheeler 
to refer informally to these collapsed stars as black holes. This evocative way of 
referring to these very strange objects soon became the standard terminology. 

In a lecture given in the memory of Karl Schwarzschild at the Astronomischen 
Gesellschaft in Hamburg, Germany, in 1986, Chandrasekhar stated the following 

Black holes are macroscopic objects with masses varying from a 
few solar masses to millions of solar masses. To the extent that 
they may be considered as stationary and isolated, to that extent, 
they are all, every single one of them, described exactly by the 
Kerr solution. This is the only instance we have of an exact de- 
scription of a macroscopic object. Macroscopic objects, as we see 
them around us, are governed by a variety of forces, derived from 
a variety of approximations to a variety of physical theories. In 
contrast, the only elements in the construction of black holes are 
our basic concepts of space and time. They are, thus, almost by 
definition, the most perfect macroscopic objects there are in the 
universe. And since the general theory of relativity provides a sin- 
gle two parameter family of solutions for their description, they 
are the simplest as well. 

The two parameter family of exact solutions referred to by Chandrasekhar is the 
family of Kerr solutions of the Einstein field equations of General Relativity, dis- 
covered by Roy Kerr in 1963, [42]. These generalize the Schwarzschild geometry to 

^The reader is referred to the excellent survey by Israel 1401 on the genesis and evolution 
of the concept of a dark star, starting from MichcU and culminating in the laws of black hole 
thermodynamics . 
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the case in which the black hole has non-zero angular momentum, and are station- 
ary and axi-symmetric as opposed to stationary and spherically symmetric. In the 
introduction to an earlier paper, written on the occasion of Einstein's centenary, 
Chandrasekhar wrote about the Kerr solution i^lDj : 

The special significance of Kerr's solution for astronomy derives 
from the theorems of Carter (1971) and Robinson (1975) which 
establish its uniqueness for an exact description of black holes that 
occur in nature. But Kerr's solution has also surpassing theoretical 
interest: it has many properties that have the aura of the miracu- 
lous about them. These properties are revealed when one considers 
the problem of the reflection and transmission of waves of different 
sorts (electromagnetic, gravitational, neutrino, and electron waves) 
by the Kerr black hole. 

Regarding the "miraculous properties" , he went on to say in the conclusion of the 
same paper: 

What, may we inquire, are these properties? In many ways, the 
most striking feature is the separability of all the standard equa- 
tions of mathematical physics in Kerr geometry. 

The goal of this paper is to survey some recent developments in the study of 
linear wave equations in Kerr geometry, and to thereby give a description of some 
of the key physical properties of black holes from a rigorous mathematical perspec- 
tive. We remark that there are also other very interesting and important results 
in mathematical relativity, like for example black hole uniqueness theorems [37], 
Hawking radiation [SH [3S] , the formation of singularities [Tl] , the positive mass 
theorem [JSl [53] and the Penrose inequalities [351 131 E] • But this paper has its main 
focus on the topics that we have worked on: the dynamics and decay properties 
of various fields in the Kerr geometry. Put differently, linear waves in the Kerr 
geometry will serve both as vehicle and guide for our voyage to black hole physics. 
As stated above by Chandrasekhar, the setting for our voyage is provided by the 
space-time geometry corresponding to an exact solution of the Einstein field equa- 
tions of gravitation, namely the Kerr solution. We shall see that the presence of 
a black hole and a space-time singularity manifests itself in remarkable and often 
unexpected ways through the behavior of waves in the space-time geometry of the 
black hole's exterior. Among the unexpected aspects in the behavior of waves, we 
shall explain that waves corresponding to localized initial data for the Dirac and 
scalar wave equations in Kerr geometry will always decay in any localized region of 
space as t tends to infinity, suggesting that the corresponding particles either "fall 
into the black hole" or "escape to infinity." Some other features are that massive 
Dirac waves decay in Kerr geometry at a rate slower than in Minkowski space, and 
the fact that spatially localized wave packets can be used to extract energy and 
angular momentum from the Kerr black hole. We shall also see that even though 
the geometric properties of the Kerr metric make it possible to solve all the known 
linear wave equations in Kerr geometry by separation of variables, it also gives rise 
to specific challenges in the analysis of these wave equations, which require the de- 
velopment of new techniques. These challenges stem from the presence of angular 
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momentum in the Kerr black hole, which in turn causes the conserved energy for 
particles and waves to not be everywhere positive. But we are anticipating what is 
to come, and before talking about black holes, we need to set the stage by recall- 
ing some of the essentials of General Relativity, and even before that, of Special 
Relativity. 

Space-time in Special Relativity is given by Minkowski space-time, that is M^, 
endowed with the standard inner product of Lorentzian signature (+, — , — , — ). Of 
particular significance are the light cones, consisting of the sets of points satisfying 
an equation of the form 

(t - t^f ~{x- x^f -{y~ y^f ~ (z - z^f = 0. (1.1) 

The light cone p.l|l describes the set of points swept out in space-time by a light 
signal emanating from the point {xq, yo, zq) at time ^o- (We always work in units in 
which the speed of light is equal to one.) The famous Michelson-Morley experiment 
showed that the speed of light is the same for all observers in uniform motion. In 
particular, the observers cannot move faster than light, and therefore their space- 
time trajectories, or world lines, lie everywhere inside the light cone. Such lines 
are called time-like. The generators of the light cone (jl.ip correspond to the world 
lines of light rays, and are referred to as null lines. The light cones thus introduce 
a causal structure in space-time, according to which every observer in uniform mo- 
tion has a future and a past, given by the interior of the upper, respectively lower 
sheet of the light cone based at its space-time position. The set of affine trans- 
formations of Minkowski space-time preserving the Lorentzian inner product (and 
in particular the light cones) forms the so-called Poincare group, the semi-direct 
product of the Lorentz group 0(1,3) by the group of space-time translations. 
All of the remarkable physical consequences of Special Relativity, such as length 
contraction, time dilation and red shift correspond to the mathematical properties 
of the Poincare group. 

While Minkowski space is a model of space-time which is very well suited to the 
study of electromagnetic phenomena such as the propagation of light, it is not suffi- 
ciently general to give a description of space-time that includes gravitation. Indeed, 
the structure of Minkowski space-time is anchored around the privileged role played 
by observers in uniform motion and therefore by linear changes of coordinates, but 
we already know from Newtonian physics that gravitation manifests itself by caus- 
ing observers in free fall to be accelerated. General Relativity emerges from Special 
Relativity by incorporating the Principle of Equivalence, which equates accelerated 
frames of reference with the presence of a gravitational field and allows in particu- 
lar for general local coordinate systems, which retain the Minkowskian character of 
space-time. In this way. General Relativity becomes a geometric theory of gravita- 
tion, in which space-time is taken to be a four-dimensional manifold M4, endowed 
with a pseudo-Riemannian metric g of Lorentzian signature (+, — , — , — ), implying 
that the tangent space at each point of space-time is isometric to Minkowski space- 
time. The world lines of test particles with non-zero rest mass acted upon only 
by gravity now become time-like geodesic curves in {M/^,g), i.e. geodesies 7 whose 
tangent vector satisfies 

< 5(7,7) = 9^j{l) Y Y , 

where according to Einstein's summation convention we sum over repeated indices. 
Likewise, light rays propagate along null geodesies with (7(7, 7) = 0. These geodesic 
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curves are the analogue of the time-Hke and null lines of Minkowski space-time. In 
a local coordinate system the geodesic equation is given by 



i 

metric and its first derivatives 



where the coefficients F j^, known as the Christojfel symbols, are built from the 



r'., = ^ .g"- + a,<?,. - a..g.,) . (1.3) 

The geodesic equation can also be written as 

V^7 = 0, (1.4) 

where V is the Levi-Civita connection, often referred to as a covariant derivative, 
defined for a general vector field X by 

d 



= (d,x^ + r^., x'^ 



where d/dx^ denotes the coordinate basis of the tangent space. The Levi-Civita 
connection is the unique metric, torsion-free connection on the Lorentzian manifold. 

In 1915 Albert Einstein discovered the relativistic field equations of gravitation, 
which account for the presence of the gravitational field through space-time curva- 
ture. To formulate these equations, we need to introduce the Riemann curvature 
tensor, which quantifies how a space-time (M4, g) deviates from being Minkowskian 
in a neighborhood of any point. The Riemann curvature tensor i?' yj. is defined by 
the relations 

valid for any vector field X. It measures the degree to which the covariant deriva- 
tives fail to commute. One can easily see from the expression of the Christoffel 
symbols ()I.3|) that the Riemann tensor is linear in the second derivatives of the 
components gij of the metric, and quadratic in the first derivatives. Its expression 
in local coordinates is quite complicated, and we will not give it here. What we 
will retain is the interpretation of the Riemann tensor as a measure of the "non- 
flatness" of a metric, by recalling the classical classical theorem of Riemannian 
geometry which says that the condition i?' j-j^ = is equivalent to {M4^,g) being 
locally isometric to Minkowski space-time. 

Gravitational forces can be understood in General Relativity as tidal forces act- 
ing on nearby point particles. Mathematically, this is made precise by the Jacobi 
equation, which shows how curvature affects the behavior of neighboring geodesies. 
Consider thus a one parameter family of time-like geodesies 7 : (s,a) 1— > 7(5, a), 
where s is arclength, and a is a parameter labeling the geodesic curves in our fam- 
ily. The vector field U := dj/ds for fixed a is the unit tangent vector along the 
corresponding geodesic, while the vector field V :— d'y/da restricted to a geodesic 
curve 7(5, ao) of our family measures the deviation between that curve and the 
neighboring geodesies. The Jacobi equation expresses the second derivative of V 
along the geodesic in terms of the Riemann curvature tensor, 

= -R\^iU'^V^U' . 
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The right side of this equation measures the tidal force for an observer moving 
afong the geodesic. 

Einstein discovered the field equations of gravitation by the requirement that one 
recovers Newtonian gravity in the non-relativistic limit (i.e. small perturbations of 
Minkowski space and particles moving slowly compared to the speed of light) , and 
taking into account the conservation laws for energy and momentum. The field 
equations read 

^ ^T^Tij. (1.5) 

The left-hand side of p.Sp is built from the Ricci tensor Rij and scalar curvature 
R, which are obtained by taking the trace of the Riemann tensor, 

Rij — i?' iij , R = R^i. 

The right-hand side of (jl.Sp is the energy-momentum tensor of the matter fields. 
The left hand-side of (|1.5p is a symmetric tensor, which is divergence- fre(3 according 
to the second Bianchi identities, a set of integrability conditions satisfied by the 
Riemann tensor. As a consequence, the energy- momentum tensor must also be 
symmetric and divergence-free; this expresses the local conservation of energy and 
momentum. The Einstein field equations need to be augmented by a set of field 
equations for the matter fields. For example, in the case of the interaction of gravity 
with electromagnetism, one obtains the source-free Einstein- Maxwell equations, 

R,j = 87:{F,kF) - ^g.^FkiF'^'), (1.6) 
ViF'^' = 0, V.Fjk + VjFm + VfeFy - , (1.7) 

where Ftj denotes the electromagnetic field tensor (and we always use the metric 
to raise or lower indices). 

The Einstein field equations (jl.Sp are a highly complicated system of non-linear 
partial differential equations, which are extremely difficult to analyze in full gener- 
ality. This holds true even in the vacuum case (Ty = 0), where they reduce to the 
vacuum Einstein equations 

R,, = 0. (1.8) 

An important analytic result on the Einstein equations was obtained in the work 
by Christodoulou and Klainerman '16j , where the nonlinear stability of Minkowski 
space is proved. This work has inspired intensive research on the analysis of the 
Einstein equations, with and without matter. Since in this article we cannot enter 
all aspects of mathematical relativity, we shall restrict attention to special solutions 
of the Einstein equation and to the analysis of linear wave equations in a given 
space-time geometry. 

The first non-flat exact solution of these equations is a one-parameter family of 
static and spherically symmetric metrics discovered by Karl Schwarzschild in 1915, 

ds^ = fl- ^] dt^ -(l- ' dr^ ~ r\de^ + sin^ e V) • (1-9) 



"^One could of course add any constant multiple of the metric tensor gij to the left-hand side 
of ljl.5|l and still obtain a tensor with zero divergence. This amounts to adding a cosmological 
constant to the Einstein equations, but this constant is usually taken to be zero in applications 
to black holes, where space-time is assumed to be asymptotically flat. 
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This metric describes a black hole of mass A/ > 0. The metric is singular at r = 
2M, but this singular behavior is just an artifact of the coordinate system being 
used. Indeed, the metric can be extended analytically across the event horizon by 
introducing the Regge- Wheeler coordinate 

u = r + 2mlog(r - 2M), 

and an advanced null coordinate 

V = t + u , 

in which the metric is given by 

ds^ = il-—jdt^-2dvdr-r^{de^+sm'^9d(p^). 

Now the locus r = 2M is a null hypersurface which corresponds to the boundary of 
the black hole and is called the event horizon. Working in these coordinates, one can 
show that future- directed time-like or null curves can only cross the hypersurface 
r = 2M from outside to inside, which explains why this hypersurface is referred 
to as an event horizon. Physically speaking, the event horizon defines the infinite 
red-shift surface as seen by a distant observer. Unlike the event horizon r = 2M, 
the locus r = is a true curvature singularity, which is not removable by any change 
of coordinates. One can for example verify that the square of the Riemann tensor 
blows up at r = 0, RijkiR"^'"' ~ M^r'^ 

According to Noether's theorem, there is a close connection between symmetries 
of space-time and conservation laws. In general relativity, this connection is made 
precise using the notion of a Killing field. A Killing vector field X is characterized 
by the so-called Killing equation 

\7iXj + S/jXi = . 

This equation implies that the Lie derivative of the metric vanishes, meaning that 
the flow of the vector field X is an isometry of space-time. In order to get the 
corresponding conservation law, one contracts the energy momentum tensor with 
the Killing field. The resulting vector field y = T'^^ Xj is divergence-free. Namely, 

div(r) - v,(r^x,) = (v,r^) x, + r-'' v,x, ^ o , 

where we applied the Killing equation and used that the energy-momentum tensor is 
symmetric and divergence- free as a consequence of the Einstein equations. Applying 
Gauss' theorem, one concludes that every boundary integral of Y vanishes, 

/ Y,i^id^i9v=0 (1.10) 
Jav 

(where V is any open subset of space-time with smooth boundary, Vj denotes a 
unit normal, and dfigy is the volume form on dV). Taking for V a set whose 
boundary consists of two space-like hypersurfaces corresponding to two time slices 
of an observer, (jl.lOp gives rise to a conservation law for a spatial integral. 

We remark that the above method is also useful in the situation when space- 
time has no obvious symmetries. In this case, one can try to find a vector field 
which is an approximate Killing field in the sense that the expression ViXj + 
V jXi vanishes up to a small error term. Then in l|1.10p one also gets error terms, 
but nevertheless one can hope to obtain useful estimates. This is the so-called 
vector field method which has been used in many papers to study the behavior of 
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waves in the Schwarzschild geometry (see for example O [171 [18] ) . An extension 
of the method to conformal symmetries goes back to Morawetz [45j. For a general 
exposition on the role of vector fields in the analysis of Euler-Lagrange systems of 
partial differential equations of hyperbolic type we refer to [15J . 

2. The Kerr Metric and the Black Hole Uniqueness Theorem 

It took nearly fifty years until an exact solution of the Einstein equations de- 
scribing the outer geometry of a rotating black hole in equilibrium was found. This 
is the Kerr metric given by 

ds^ = ^{dt-asm^dd^)^-U (^^ + d^'^^-^^(^adt-{r^ + a^)d^y, (2.1) 

where 

U = + a^cos^'d, A = - 2Mr + a^ (2.2) 
and the coordinates {t, r, ip) are in the range 



-oo<t<oo, M + - a2 < r < oo, < < tt, < < 27r. (2.3) 

Here the parameters M and a describe the mass and the angular momentum per 
unit mass of the black hole. It is easily verified that in the case a = 0, one recovers 
the Schwarzschild metric (|1.9p . 

For this metric to describe a black hole we need to assume that > , giving 
a bound for the angular momentum relative to the mass. In this so-called non- 
extreme case, it can be verified by an argument similar to the one used for the 
Schwarzschild solution that the null hypcrsurfacc 



r = ri := M + \J - (2.4) 

defines the event horizon, the boundary of the black hole. Here we shall only con- 
sider the region r > r\ outside the event horizon. The coefficients of the Kerr metric 
are independent of t and (/?, showing that the space-time geometry is stationary and 
axi-symmetric. One of the key features of the Kerr geometry is the existence of 
an ergosphere, that is a region which lies outside the event horizon, in which the 
vector field ^ is space-like. In order to determine the ergosphere, we consider the 
norm of the vector field — 

, A-a^sin^i? r"^ - 2Mr + cos^ 'd . , 

9^] = 9tt = jj = jj . (2.5) 

This shows that ^ is space-like in the open region of space-time where 

- 2Mr + cos^ z9 < , (2.6) 

which defines the ergosphere. It is a bounded region of space-time outside the event 
horizon, which intersects the event horizon at the poles = 0, tt. The implications 
of the ergosphere for the analysis of wave equations cannot be overstated because 
the ergosphere corresponds to a region in which the conserved energy associated to 
a given field will fail to he positive (see for example (j5.2[) and (|5.3p below) . 

The physical significance of the Kerr metric manifests itself through the black 
hole uniqueness theorems of Israel [38] , Carter [8] and Robinson |47l [37] . Physically, 
these theorems indicate that the stationary end-state of a rotating black hole should 
be given precisely by the non-extreme Kerr geometry. From the mathematical point 
of view, these theorems are uniqueness results for a class of boundary value problems 
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for the vacuum Einstein equations, and their proof requires additional mathemat- 
ical assumptions. As boundary conditions one assumes the asymptotic flatness of 
space-time in the form of weak asymptotic simplicity (see [52' page 282]), plus 
the existence of an event horizon with the natural spherical topology. Furthermore, 
one assumes axi-symmetry and pseudo-stationarity (pseudo-stationarity means that 
there is a Killing field which is assumed to be timelike only in the asymptotic end). 
Finally, as technical assumptions which are also motivated from physics, one needs 
to impose a causality condition as well as time orientability, and one needs to fix 
the topology of space-time to be x S*^. (The review article by Carter [5^ and the 
books [33l [37] specify and discuss all the conditions in detail.) The theorem is as 
follows: 

Theorem 2.1. Under the above assumptions, every solution of the vacuum Ein- 
stein equations admits a global chart {t, r, "iJ, ip) in which the metric is the Kerr 
solution ()2.ip . 



3. Linear Wave Equations in the Kerr Geometry and their Separation 



In this section we give an overview of the physically relevant linear wave equations 
and present their separability properties in the Kerr geometry. For simplicity, we 
begin with the scalar wave equation, which in a general space-time reads 

d 



1 



-g dx^ 



0, 



(3.1) 



where g denotes the determinant of the metric gij . For the Kerr metric, this becomes 



— A — 

dr dr 

d 



1 

A 



d_ 

dip 



dcos-d 



sin ?9 



d 



1 



d cos -d siv? § 



asm^^— + — 
at dip 



$ = 0. 



(3.2) 



We denote the square bracket in this equation by □, (although it is actually a 
scalar function times the wave operator in (|3.ip ). We now explain the separation 
of variables as discovered by Carter [3] . Due to the stationarity and axi-symmetry, 
we can separate the t- and (^-dependence with the usual plane-wave ansatz 



$(t,r,i?,(^) 



-iut—ikip 



(3.3) 



where w is a quantum number which could be real or complex and which corresponds 
to the "energy" , and k G Z is a quantum number corresponding to the projection of 
angular momentum onto the axis of symmetry of the black hole. Substituting (|3.3p 
into (j3.2p . we see that the wave operator splits into the sum of radial and angular 
parts, 

= {n^.k+A^,k)^, (3.4) 
where Tiuj,k and Auj^k are given by 

1 



T^LU.k 

■Au.k 



—A— 

dr dr 

d 

dcosd 



-{{r^ +a^)LU + aky, 



sin^ •& 



d 



1 



d cos sin^ i? 



{au sin^?9 + fc)' 



(3.5) 
(3.6) 
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We can therefore separate the variables r and •& with the muhiphcative ansatz 

0(r,^) - i?(r)e(z9), (3.7) 

to obtain for given uj and k the system of ordinary differential equations 

7^c.,fei?A = -Ai?A, A^,kQx = AOa. (3.8) 

The separation constant A is an eigenvalue of the angular operator A^^k and can 
thus be thought of as an angular quantum number. In the spherically symmetric 
case (i.e. a = 0), A goes over to the usual eigenvalues A = Z(/ + 1) of the total angular 
momentum. We point out that the last separation (13. 7p is not obvious because it 
does not correspond to an underlying space-time symmetry. 

For the Dirac equation describing a spin i field, the situation is more compli- 
cated because it is a system of differential equations and involves a rest mass m. 
The separability of the Dirac equation in the Kerr metric was first established by 
Chandrasekhar in 1976, using an ingenious new method. The Dirac equation 
for a particle of mass m reads 

(?7^ Vj - m)-^ = 0, (3.9) 

where the Dirac matrices 7^ are related to the metric by the anti-commutation 
relations 

and V is a connection on the spinors. We write the Dirac equation in a Newman- 
Penrose null frame IllJ (/,n,m,m), i.e. in a frame in which the metric takes the 
form 

Qij = liiiij -\- Uilj — mifhj — fhiTrij . 

Choosing a suitable frame and after an appropriate linear transformation of the 
Dirac spinor, 

where S — S(r, 1!)) is a diagonal matrix, the Dirac equation becomes 

(7^ + y^)^ 0, (3.10) 

where TZ and A are certain first-order matrix differential operators. The separa- 
tion of variables is achieved by assuming that each component of the transformed 
spinor 4* is a product of functions of one variable, according to a specific separation 
pattern, 

X+(r)y_(z9) 

By substituting (|3.1ip into (|3.10p . we obtain the separated matrix equations 



keZ. (3.11) 



7^* = A*, y^* = -A*, (3.12) 

where A is a separation constant, under which the transformed Dirac equation (|3.10[) 
decouples into a system of ODEs. 

The remaining field equations of physical interest are those describing electro- 
magnetic and linearized gravitational waves. It is a remarkable and very useful fact 
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that these waves are all governed by a single second-order equation, the so-called 
Teukolsky master equation It reads 



-A- - ^ (.^ + «^)^ + - - - 4.(. + .acos^)- 



■ sm^ 6*- H ^ a -svcl ^ 



9 cos d d cos sin^ 6 \ dt 



(3.13) 



= 0. 




The parameter s denotes the spin. We explain how this equation comes about, 
beginning with the case s = ±1 of a Maxwell field. The source- free Maxwell's 
equations are given by 

V/F'^' = 0, y^Fjk + VjFfe, + VfcFy = , (3.14) 

where Fij is the electromagnetic field tensor. Choosing the same Newman-Penrose 
null frame as in the separation of the Dirac equation, we can combine the compo- 
nents of the field tensor into the three complex functions 

where p ^ ~{r ~ ia cos S)^^ . Then Maxwell's equations (|3.14p can be written as 

*i 

where P is a first-order matrix differential operator. The key point is that, multiply- 
ing this equation from the left by a suitable first-order matrix differential operator, 
one obtains decoupled second order equations for ^E"! and ^'-i, which are precisely 
the Teukolsky equation for s = ±1. It is important to note that this process of 
cross-differentiation and elimination makes use of some key commutation identities 
between certain covariant derivative operators in the Kerr metric. These identities 
hold as a consequence of the special algebraic structure of the Riemann tensor in the 
Kerr geometry; namely, it is of type D in the Petrov-Penrose classification, mean- 
ing that the Weyl tensor has two repeated principal null directions. If 5*1 or 
are known, then the remaining components are readily obtained by differentiation 
using the so-called Teukolsky-Starobinsky identities [TT] . 

Finally, for s = ±2, the Teukolsky equation (|3.13p is obtained similarly by cross- 
differentiation and elimination on the systems of first-order equations obtained by 
linearizing the Bianchi identities 

^iRjklm + ^jRkilm + ^ kRijlm — 0, 

around the Kerr metric. If we denote by V'o^'' and ip'^'^ the first-order perturbations 
of the Newman-Penrose curvature components t/^o and ip^ defined by 

where Cijki denotes the Weyl conformal curvature tensor, then the unknowns in 
the Teukolsky equations with s = ±2 are given by 

We remark that for s — 0, the Teukolsky equation (|3.13p reduces to the scalar 
wave equation p.ip . The case s = 1/2 corresponds to the massless Dirac equation. 
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We also point out that the reduction to a second order scalar equation does not work 
for the massive Dirac equation. Just like the scalar wave equation, the Teukolsky 
equation can be reduced to ordinary differential equations by the separation ansatz 

- e-"^'-'^"^ Rs{r) e,(i?) . 

These are the separability properties that were referred to by Chandrasekhar [TU] 
in the excerpt quoted in the Introduction as "having the aura of the miraculous." 

Clearly, the above linear wave equations can also be analyzed in the special 
case a = 0, where the Kerr geometry simplifies to the spherically symmetric and 
static Schwarzschild geometry (jl.9p . Since in this paper we shall concentrate on 
rotating black holes, we merely mention a few important papers on the analysis 
of waves in the Schwarzschild geometry. In the fundametal paper [41j it is proven 
that the solutions of the scalar wave equation are pointwise bounded, uniformly in 
time. The corresponding rates of decay are derived in [12] and [121 HI]- In [5] the 
boundedness of solutions to Maxwell's equations in the Schwarzschild geometry is 
proved. Finally, |20| considers the scalar wave equations in the Kerr metric in the 
case a <ti M of small angular momentum, where the metric is close to Schwarzschild. 

4. Dynamics of Dirac Waves 

To get a more detailed picture of the properties of a rotating black hole, it is very 
helpful to consider dynamical situations where physical objects are moving in the 
Kerr geometry. The simplest situation is to consider point particles, whose motion 
is described by the geodesic equation (II. 2p . As is worked out in detail in [TT| , there 
are stable orbits where a point particle "rotates around" the black hole. Apart 
from the circular orbits, there are elliptic orbits which are not quite closed due to 
the perihelion shift. Furthermore, there are orbits where the particle falls into the 
black hole or escapes to infinity. Considering quantum mechanical particles makes 
the situation more interesting, because due to the Heisenberg Uncertainty Principle, 
these phenomena could happen simultaneously with certain probabilities. To clarify 
the picture, in this section we consider the dynamics of Dirac waves in the Kerr 
geometry. We shall first see how the Cauchy problem for the Dirac equation in Kerr 
geometry can be solved by means of an integral representation for the propagator, 
realized as a superposition of "modes" arising from the separation of variables. We 
then analyze this integral representation to show that in the presence of a black hole, 
a quantum mechanical Dirac particle cannot remain in a localized region of space 
for arbitrarily large times, in other words that the quantum mechanical particle 
corresponding to the Dirac wave function will either eventually "fall into the black 
hole" or "escape to infinity" . We also discuss the decay rates. Our analysis is 
of course limited by the fact that we look at a test Dirac field in Kerr geometry 
as opposed the fully coupled axisymmetric Einstein-Dirac equations, but the non- 
existence results proved in [55] for black hole solutions of the spherically symmetric 
Einstein-Dirac-Maxwell equations suggest that the same conclusions should hold in 
the fully coupled axisymmetric casejj. 

The first step is to bring the Dirac equation ()3.10p into Hamiltonian form, in a 
way that is compatible with the separation of variables. This is not a trivial step 

major obstacle in the analysis of the fully coupled axisymmetric Einstein-Dirac-Maxwell 
equations by similar methods is that a complete separation of variables into ordinary differential 
equations seems most unlikely, and one has to deal with the full system of PDEs. 
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because the coefRcients of the ft-term depend on both r and ^. The Dirac equation 
becomes 

d 

f — * = i7 * , (4.1) 
dt ^ ^ 

where i7 is a first-order matrix differential operator, ft is important for the analysis 

of this evolution equation to introduce a positive scalar product <.|.>, with respect 

to which H is symmetric. To this end, we consider the Dirac current 5'7-'^', which, 

as a consequence of the Dirac equation, it is divergence free. Furthermore, it is 

timelike and future-directed, so that its inner product \I>7°v[/ with the unit normal 

to the space-like hypersurface t = const is positive, even inside the ergosphere. We 

introduce |$> by integrating the corresponding bilinear form over the 

hypersurface t — const. Then the Gauss divergence theorem yields that <5' | <i>> 

is time independent. Hence, using (|4.ip . 

= 4<*(*)l*(^)> = i(<H^{t)\<S>(t)>-<^{t)\H^(t)>), 
dt 

and thus H is indeed a symmetric operator on the corresponding Hilbert space Ti.. 

After constructing a self-adjoint extension, we can use the functional calculus 
and the spectral theorem in Hilbert spaces to obtain 

^{t) ^ e-^^"^o ^ [ e-'"* dE^ . 

Ja(H) 

Using the separation of variables, we can analyze the spectrum and the spectral 
measure. More precisely, as a technical tool for analyzing the spectral measure, we 
first set up the problem in a finite radial box, with Dirichlet boundary conditions 
on the walls of the box that make the Dirac operator essentially self-adjoint. The 
Hamiltonian for this problem has a purely discrete spectrum and one thus obtains a 
spectral representation of the propagator as an infinite sum of discrete eigenstates, 
expressible as products of eigenfunctions of the radial and angular operators arising 
from the separation of variables. The next step is to take limits as the walls of the 
box tend to the event horizon and to infinity. In these limits, the sum over the 
discrete eigenvalues goes over to an integral representation for the solution of the 
Cauchy problem. We thus obtain the following integral representation. 

Theorem 4.1. Consider the Cauchy problem for the Dirac equation in the Kerr 
geometry, 

{i-f^Vj-m)^{t,x) = 0, *(0,a;) = *o(a;), (4.2) 
for initial data *o G Cq' {{n, oo) x S'^)^. Then 

-1 r + oo 2 

'^it,x) ^ - J2 J2 t^r*a'^"(a;)<«'b"" I «'o>rfw, (4.3) 

k,neZ''~°° a, 6=1 

where the sums and the integrals converge in the Ti-norm. 

In the above integral representation, the integer k is the quantum number corre- 
sponding to the projection of the angular momentum on the axis of symmetry of 
the black hole, and n is the generalized total angular momentum quantum number 
arising from the separation of variables, corresponding to the discrete spectrum of 
the angular part of the Dirac operator in the Kerr geometry. The are solu- 

tions of the Dirac equation which arise from the separation of variables and satisfy 
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asymptotic boundary conditions at infinity and near the event horizon. More pre- 
cisely, for > TO the ^|;^„ behave near the event horizon hke incoming spherical 
waves for a = 1 and outgoing spherical waves for a = 2, whereas in the case \uj\ < to, 
the vj/^"" is a linear combination of both incoming and outgoing spherical waves 
near the event horizon, such that the solution has exponential decay at infinity. 
The coefficients t^^" are given explicitly in terms of the reflection and transmission 
coefficients. 

Applying the Riemann-Lebesgue lemma to the above integral representation 
gives the following local decay result. 

Theorem 4.2. Consider the Cauchy problem ()4.2|) for the Dirac equation in the 
Kerr geometry, where the initial data ^'o is in L'^((ri,oo) x S'^)* and in near 
the horizon, i.e. \'^o{x)\ < c for x g (r\,r\ + e) x S*^. Let S > be given, let 
R > ri + 6 and consider the compact space-like hyper-surface Kg.R of £ given by 

Ks.R {(t, r,e,(j))\ri+5<r<R,t^ con.st.}. (4.4) 

Then the probability for the Dirac particle to be inside Kg.R tends to zero as t oo, 
that is 

lim / (^-f°^){t,x)dn = 0, (4.5) 

where dfj, denotes the induced volume element on Kg R. 

Theorem 14.21 implies that the Dirac spinor 5* decays to zero in L"^^, or equiva- 
lently that the Dirac particle must eventually either disappear into the black hole, 
or escape to infinity. In order to get a more detailed physical picture, one would 
like to determine the probability of these outcomes in terms of the Cauchy data. 
Likewise, one would like to determine the rates of decay of the Dirac spinor in Lf^^, 
as t tends to infinity. Both of these questions were addressed in [2^, under the 
additional assumption that in the integral representation (|4.3p only a finite number 
of angular momentum modes are present, 

|fc| < fco , |"| < no . (4.6) 

We now state the main results, beginning with the decay rates. 

Theorem 4.3. Consider the Cauchy problem as in Theorem \4-.2\ with initial data 
normalized 6y <^'o | 5'o> = 1- Suppose that (|4.6p holds. Then we have: 

i) If for any k and n, 

limsup|<«'^'^" I *o>| 7^ 0, or liminf |<*^'^" | ^'o>| 7^ 0, (4.7) 

then, as t ^ 00, 

\^{x, t)\ = ct-'^^^ + 0(f-5/6-^)^ (4 g) 

where c = c{x) =/= 0, and e < 1/30. 

a) If for all k, n and a — 1,2, <vlijj'^" | vI/q> = Q for all lo is a neighborhood of ztm, 
then for any fixed x, ^(x, t) decays rapidly in t. 

Note that the decay rate obtained in Theorem 14.31 is slower than the rate of decay 
of t^'^/^ one obtains for the solutions of the Dirac equation in Minkowski space, |26) . 
At first glance, this result seems to be somewhat counter-intuitive. Indeed, because 
the Kerr metric is asymptotically flat, the Dirac spinor should behave near infinity 
like a solution of the Dirac equation in Minkowski space, where the rate of decay 
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is in t^^^"^. On the other hand, near the event horizon, the Dirac particle should 
behave like a massless particle, and its wave function should have rapid decay. One 
would therefore expect the rate of decay in the Kerr metric to be at least as fast as 
^-3/2_ ^jj^jg naive picture is in fact incorrect. Indeed, as shown in [26j, the slower 
rate of decay can be understood by the fact that in the Kerr geometry, the energy 
spectrum of the initial data for the Dirac equation oscillates more and more as m 
approaches the rest mass m of the Dirac particle (essentially as sin((TO — 
as opposed to {uj — m)^/^ in Minkowski space). Upon taking the Fourier transform 
in time, these oscillatory contributions lead to a slower decay rate. Another way 
of understanding this effect is that the backscattering of the outgoing wave near 
infinity slows down the decay. 

We now turn to the probability estimates. The probability for the particle to 
escape to infinity is given by 



p = lim / ■9y'$(t,x) i/j dfj,, (4.9) 

Jr>R 

where R > ri. In [26 it is shown that the probability p is independent of R and 
can be expressed in terms of the quantities appearing in the integral representation 
gS) by 



Analyzing this expression one gets the following result. 

Theorem 4.4. Consider the Cauchy problem as in Theorem \4-.2\ with initial data 
normalized by <^'o | 5'o> = 1- Then the following hold: 

(i) // the outgoing initial energy distribution satisfies <\I'|'^" | 5'o> ^ for 
some uj such that \uj\ > m, then p > Q. 

(ii) // the initial energy distribution satisfies for a = 1 or a — 2, <vl/^'^" | vI/q> ^ 
for some lo such that \lu\ > m, then p < I. 

(iii) If the initial energy distribution is supported in the interval [—m,m], then 
p = 0. 

(iv) If for any k and n, 

limsup|<*^'^" |*o>| 7^0 or liminf |<*^'^" | ^'o>| 7^ 0, (4.10) 
then <p < 1. 

(v) We have p ^ 1 if and only if for all k and n the following conditions hold: 

<^rfc-" I = 0, |c^|<m, (4.11) 

<^J'^» I *o> = -2t5;'^"<*^'^" I *o> ^f |w| > m. (4.12) 

We point out that in the above theorems we always considered solutions sup- 
ported outside the event horizon. One can also analyze the Dirac equation across 
the event horizon in the weak sense, and show that there are no time-periodic weak 
solutions which are in away from the event horizon [23l |24] . 

We finally mention a few other rigorous results on the dynamics of quantum 
mechanical waves in the Kerr geometry which we we cannot describe in this article. 
First, the scattering theory for Dirac particle has been developed in [36l [211 E]- 
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A remarkable result on the Klein-Gordon equation was obtained in 34J. The pa- 
per is devoted to the Hawking effect, which leads to the creation of fermions 
by a rotating black hole. 

5. An Integral Representation and Decay for Scalar Waves 

For the scalar wave equation, we derive an integral representation which is similar 
to that obtained for the Dirac operator in Theorem l4.H but the proof is significantly 
more difficult and requires new techniques. The integral representation will imply 
decay as for the Dirac equation, but it will in addition lead to the phenomenon 
of superradiance, which will be explained in Section [SI In contrast to the Dirac 
operator, the scalar wave operator in Kerr geometry does not admit a conserved 
quantity which is positive everywhere outside the event horizon, making it impos- 
sible to apply the spectral theory of self-adjoint operators in Hilbert space. First 
of all, the energy E[^] is not positive. Namely, using the invariance of the scalar 
wave Lagrangian 

= |V$P (5.1) 
under time translations, we obtain by Noether's theorem the expression 

£:[$] = / dr d(cosz9) / -^£, (5.2) 



2tt 



where £ is the energy density 



(r2 + a2)2 



,2.q|Q _^|2 / 1 



+ sin^^9|9cos^*r + X l^^^l - (5-3) 

One sees that all the terms in the energy density are positive, except for the co- 
efficient of |i9i^$p, which is positive if and only if — 2Mr + cos^ d > 0, that 
is, precisely outside the ergosphere. The difficulty is that the conservation of E[^] 
does not give a Sobolev estimate for <I>. In particular, energy conservation does 
not rule out the situation where $ blows up in time, in such a way that the energy 
density tends to minus infinity inside the ergosphere and to plus infinity outside the 
ergosphere. Moreover, in contrast to the Dirac equation, also the conserved charge 
is non-positive. In |27| it is shown that no other first order or higher order positive 
conserved energy exists for the scalar wave equation in the Kerr geometry. 

We begin by stating a theorem which gives an integral spectral representation 
for the solution of the Cauchy problem for the scalar wave equation. This theorem 
can be thought of as the analogue of Theorem 14. II for the scalar wave equation. 

Theorem 5.1. Given initial data 'i'o e C^(M x 5^)^, the solution of the Cauchy 
problem for the scalar wave equation (|3.2p can be represented as 

1 r ^ /I 

(5.4) 



where 



— uj — Lo^^ and = 5 9 ■ (5-5) 



Here the sums and the integrals converge in Lj^ 
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We now describe the various terms that appear in the above integral representation 
in some detail, since they will be used in our discussion of superradiance in Section 
ini Recall from the discussion of the separation of variables for the scalar wave 
equation that with the separation ansatz (j3.3[) , the scalar wave operator splits into 
the sum of a radial operator (j3.5[) and an angular operator (j3.6l) , with separated or- 
dinary differential equations (|3.8p . The angular operator A^j^k has a purely discrete 
spectrum of non-degenerate eigenvalues < Ai < A2 < • • ■ . The corresponding 
eigenfunctions Of^''' are the spheroidal wave functions. In order to bring the radial 
equation into a convenient form, we introduce a new radial function (j){r) by 



2 , Viu)]^iu) - 0. (5.8) 



<j){r) = \Jr'^ + a2 R{r) , (5.6) 
and define the Regge- Wheeler variable m G M by 

du + ,^ 

mapping the event horizon to u = —00. The radial equation then takes the form of 
a Schrodinger equation, 

du 

For fixed separation constants fc > 0, w and An, the potential V{u) has the asymp- 
totics 

lim V{u) ^ , lim V{u) ^ ~uj^ . 

u — * — CKj u — *oo 

Thus there are fundamental solutions and </) of (|5.8p which behave asymptotically 
like plane waves at the event horizon and near infinity, respectively, 

(f>{u) ^ e'"" as u ^ -00, 0(u) ~ g-*'^" as u — > 00 

(these functions are the so-called Jost solutions, see [28] for details). As the corre- 
sponding time dependent solutions behave in time like the plane wave e~'"*, these 
solutions have a physical interpretation as incoming and outgoing waves in the as- 
ymptotic regions. The functions (p and (f> form a fundamental system, and thus we 
can represent (j) as 

(j){u) = A^{u) + B^{u) . (5.9) 

The coefficients A and B are the reflection and transmission coefficients. The quan- 
tities tab in the integral representation (|5.4p are explicit functions of A and B. Fi- 
nally the functions '^kuni^^ a = 1, 2 are the solutions of the wave equation (13. 2|) . 
with fixed quantum numbers k,uj,n, corresponding to the real- valued fundamental 
solutions of the radial equation given by 

(/)^=Re(^, 02=:Im(^. (5.10) 

Our next result is a decay theorem for scalar waves in Kerr geometry analogous 
to Theorem 14.21 for the Dirac operator. 

Theorem 5.2. Consider the Cauchy problem for the wave equation in the Kerr 
geometry for smooth initial data which is compactly supported outside the event 
horizon and has fixed angular momentum in the direction of the rotation axis of the 
black hole, i.e. for some k G 'Z, 

(a>o,9t<i>o) = e-''=^($o,9t$o)(r,i?) e {{n , 00) x . 
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Then the solution decays in L'^^{{ri,oo) x S'^)^ as t —> oo. 

Theorem 15.21 can be interpreted as a linear stability result for the Kerr black hole 
under perturbations by massless scalar fields for a finite number of azimuthal an- 
gular momentum modes. It would be desirable to analyze the convergence of the 
sum over all /c-modes. For this, one must obtain estimates which are uniform in k. 
Such estimates have not worked out, except in the case a <C m, a perturbation 
of the Schwarzschild space-time; see [20], where the authors obtain a boundedness 
result using the vector field method. No theorems are known on rates of decay 
or probability estimates for scalar fields that would be analogous to Theorems 14.31 
and 14.41 for Dirac fields. Proving such theorems in Kerr geometry appears to be a 
significant challenge. 

We conclude this section by sketching the proof of Theorem 1 5. II As for the Dirac 
equation, we reformulate the wave equation in Hamiltonian form. Letting 

* = ( ) ' ^'-''^ 

the wave equation (|3.2p takes the form 

idt'i! = H^, (5.12) 

where H is the Hamiltonian 

and where a and /3 are certain differential operators. The general strategy for 
constructing an integral representation for the propagator is the same as in the 
case of the Dirac operator, and requires as a first step that the problem be set 
up in a finite radial box, with Dirichlet boundary conditions on the walls of the 
box. Unfortunately, the lack of a positive conserved energy makes it impossible 
to associate a positive definite invariant inner product to the time evolution of 
the Hamiltonian, and one has to consider the Hamiltonian H as acting on a func- 
tion space endowed with an inner product of indefinite signature. In the case of 
a finite box, the inner product space is a Pontrjagin space [61, that is a complex 
vector space /C endowed with a non-degenerate inner product <. , .> and an or- 
thogonal direct sum decomposition K. = /C+ ® /C_, such that (/C+,<. , .>) and 
(/C_, — <. , .>) are both Hilbert spaces, with /C„ being finite-dimensional. Classical 
results of Pontrjagin [6] imply that any self-adjoint operator A on a Pontrjagin 
space has a spectral decomposition, which is similar to the one given by the spec- 
tral theorem in Hilbert spaces, except that there is in general an additional finite 
point spectrum in the complex plane (which is symmetric about the real axis). We 
now consider the vector space Vr^.m — {H^'^ ® L^)([ri,rfl;] x S*^) with Dirichlet 
boundary conditions ^i(ri) = = ^'i(rfl). We endow this vector space with the 
inner product associated to the energy ()5.2|) . It can be shown [27] that for every 
fR > ri there is a countable set E C {ri,r]i) such that for all g (?'i,7'_r) \ E, 
the inner product space Vrj^^m is a Pontrjagin space with the same topology as 
(i?^'^®L^)([ri, Tfl] X 5^), and that the Hamiltonian on the Pontrjagin space 'PrL,rii 
with domain V = C°°{[rL,rii] x 5^)^ C 'PrL,rR is essentially self-adjoint. If the size 
of the box increases, the number of complex conjugate pairs becomes larger, and 
in the infinite volume limit all these spectral points move onto the real axis. This 
is made precise using contour methods and estimates for the resolvent. Whiting's 
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Figure 1. Superradiant scattering. 



mode stability result [53], which guarantees the absence of unstable exponentially 
growing modes for the separated wave equation, is crucial to our proof. 

6. A Rigorous Treatment of Superradiance for Scalar Waves 

One of the most fascinating aspects of the classical physics of black holes is given 
by the so-called Penrose process [46^ , which shows that one can extract energy and 
angular momentum from the Kerr black hole so as to lower the angular momentum 
to zero by suitably exploiting the effect of the ergosphere on the dynamics of point 
particles. The basic idea of the Penrose process is as follows (see [52] for more details 
and [51] for more realistic scenarios involving collisions of charged particles). First 
recall that the conserved energy of a point particle of momentum p is given by {p, dt) , 
and that this inner product clearly need not be positive inside the ergosphere, 
where the vector field dt is space-like. Now consider a rocket which flies into the 
ergosphere, where it splits into two objects whose energies have opposite signs. 
By a suitable choice of the energy and momenta, one can arrange that the object 
of negative energy crosses the event horizon and reduces the energy and angular 
momentum of the black hole, whereas the other object escapes to infinity, carrying 
(due to energy conservation) more energy than the original rocket. In this way, 
one can extract energy from the black hole, at the expense of reducing its angular 
momentum. Christodoulou 13J showed that the infinitesimal changes of mass 5M 
and angular momentum 5{aM) of the black hole satisfy the inequalities 

5[aM) < d±^SM < 0, (6.1) 
a 

and as a consequence he showed that it is not possible to reduce the mass of the 
black hole via the Penrose process below the so-called irreducible mass 

Mf„ i (Af2 + ^M^-{aMY) . (6.2) 

Superradiance is the wave analogue of the Penrose process. One considers a wave 
entering the black hole. One part of the wave enters the event horizon, whereas the 
other part is scattered at the black hole and gives rise to an outgoing wave (see Fig- 
ure [1]). By arranging the energy of the infalling wave to be negative, one can again 
extract energy from the black hole. This effect was first studied by Zel'dovich [55] 
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and Starobinsky |49j on the level of modes, i.e. by considering the transmission and 
reflection coefficients in (|5.9p associated to the Schrodinger equation ()5.8|) . In this 
analysis, the quantities uP'\A\^ and w^jiJp have the interpretations as the energy 
flux of the incoming and outgoing waves, respectively. Thus the relative energy 
gain [H is given by 

^ = ^ • (6.3) 

Computing the Wronskians of and near the event horizon and near infinity 
gives the relation 

- - . (6.4) 

ijj 

If the right side of (|6.4p is positive, the outgoing flux is smaller than the incoming 
flux, and this corresponds to ordinary scattering. However, if the right side of (|6.4p 
is negative, then the outgoing flux is larger than the incoming flux, and accord- 
ing to (16. 3|) we gain energy. This is termed superradiant scattering. Using (j5.5|) . 
superradiant scattering appears precisely when co is in the range 

< < Ic^ol . (6.5) 

Starobinsky [49] computed 91 and found a relative gain of energy of about 5% for 
k = 1 and less than 1% for k > 2. Teukolsky and Press [50] made a similar mode 
analysis for higher spin and found numerically an energy gain of at most 4.4% for 
Maxwell (s — 1) and up to 138% for gravitational waves (s = 2). For more recent 
developments on the mode analysis of superradiance and related physical effects we 
refer to [7j. 

Unfortunately, the mode analysis does not give information on the dynamics. 
Thus for a rigorous treatment of energy extraction one needs to consider the time- 
dependent situation. This was done numerically in [T] for wave packet initial 
data. The main result of [29] treats the time-dependent situation rigorously for 
the Cauchy problem. As initial data we take wave packets of the form 



\ I iuju ( 1 

+ Cout e 



where L is a large parameter, and 77^ is a smooth cutoff function of the form 

1 (u-l? 
r]L{u) = —^ri 



(6.6) 



VI'V L 

with rj G C^(R+). Here 6n,(i(f9) is an eigenfunction of the angular operator A. 
The energy radiated to inflnity is defined by 

£:out = lim <«'(i),X(2n.oo)(0*(0> , (6.7) 

where x is the characteristic function. 

Theorem 6.1. For any R > ri and S > there is initial data <E C^{{R, 00) x 
5^)^ of the form (|6.6p such that the limit in (|6.7p exists and 



<*o,^'o> 



< 6 



with d\ as in (|6.3p . 
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For the proof we consider the integral representation of Theorem 15.11 for the wave 
packet initial data ()6.6p . The crucial analytical ingredient in the proof is the time- 
independent energy estimate for the outgoing wave as derived in ^31] . 

It should be stressed that we allow the initial data to be supported arbitrarily 
far away from the event horizon. This is important in order to avoid artificial initial 
data which would not correspond to an energy extraction mechanism. For example, 
if one allows the support of the initial data to intersect the ergosphere, one could 
take initial data with zero total energy, in which case the quotient i?out/ <^'o7 ^o> 
could be made arbitrarily large. 



7. The Stability Problem for Kerr Black Holes 

The remaining challenge is to prove the linear stability of the Kerr black hole 
under electromagnetic and gravitational perturbations. As Frolov and Novikov put 
it [321 S. 143], this is "one of the few truly outstanding problems that remain in 
the field of black hole perturbations." Since the linear perturbations of the Kerr 
metric are described by the solutions of the Teukolsky equation (|3.13p . proving 
linear stability amounts to showing that the solutions of the Cauchy problem for 
the Teukolsky equation for compactly supported Cauchy data decay in as t 
tends to infinity. In other words, the task is to prove the analogue of Theorem 15.21 
for the solutions of the Teukolsky equation (j3.13p in the cases s = 1 and s = 2. The 
analysis of the Teukolsky equation for higher spin appears to be considerably more 
difficult than that of the scalar wave equation. This is because, in contrast to the 
scalar wave equation (j5.ip . the Teukolsky equation has no variational formulation, 
and thus there is no simple method for deriving the conserved quantities. Indeed, 
the expressions for the physical energy and charge are very complicated and involve 
higher derivatives of the field ^'s. The only result so far on the Teukolsky equation 
in Kerr is the deep and important mode stability theorem proved by Whiting in '54] . 

More recently, a rigorous stability result was proved for the Teukolsky equation 
in the simpler setting of the Schwarzschild geometry |30j for electromagnetic and 
gravitational wave perturbations, which we now describe. The evolution of a mass- 
less wave of spin s in the Schwarzschild geometry is described by the Teukolsky 
equation (set a = in (|3.2p ') 

drAdr - ^ {r^dt - (r - M)s)^ - Asrdt 



-|-9cos,5 sin^ d dcos§ + . \ ^ {du, + is cos t?)^ 
sm if 



$(t,r,^,^) = 0, (7.1) 



where A = 2Mr. We consider (m|) with initial data 

$lt=o = $0 , 9t$|t=o = *i • (7.2) 
Then the following theorem holds (see [30]). 

Theorem 7.1. For spin 5 = 1 or s ^ 2, the solution of the Cauchy problem (17. ip . 
(Ol) /or ($o,*i) e C(f ((2A/,c») X S'^f decays m L^^{{2M,oo) x 5^) ast^oo. 



The proof has some novel features, and we shall discuss some of them. 
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After separating the time and angular dependence similar to (|3.3p , (|3.8p (with k £ 
Z, A G M and u; g C) and replacing r by the Regge- Wheeler variable u, the Teukol- 
sky equation can be reduced to a one-dimensional Schrodinger equation 

-^Hu) + Viu)<jyiu) = 

with potential V given by 

(r - M)2 s2 52^ A 

^^4 ^ + ^ ' 

where A is an eigenvalue of the angular operator. Due to the presence of the dt- 
term in (|7.ip . V is complex even for real lu. Thus most standard techniques are 
unavailable. 

First one easily checks that 

lim V(u) ~ — i LU 

^ ' \ AM J 

Then writing the time-dependent equations in Hamiltonian form, 

= , ^ = ($,at$) , 

with H a matrix differential operator, we can show that if lu is outside the strip 
< luiLu < then lu lies in the resolvent set of H, and the resolvent R^^ 

is holomorphic there. By means of suitable resolvent estimates, we show that 
any 'if G has the integral representation 

"ifiu) = --^ lim / {R^-9){u) diu , 

where Cn is the contour shown in Figure[21 which has two connected components Ci 
and C2. Using contour deformation techniques, we can prove the following non- 
standard integral representation theorem for the solution. 

Theorem 7.2. For spin s — \ or s — 2, the solution of the Cauchy problem for 
the Teukolsky equation (|7.ip with initial data ^E'o = {^,dt^)\t=o G C^{M.)'^ has the 



Viu) 



2{r-M) 4A 
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representation 

where is the limit of the integral kernel of the resolvemt from the lower half 
plane, 

^wo(u,w) := lim Ruj{u,v) . 
Both integrands in the above integral representation are in . 

Decay is an immediate consequence of this theorem: For t —^ oo, the first integral 
tends to zero by the Riemann-Lebesgue lemma, while the second integral tends 
exponentially to zero. 
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